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Abstract — Numerical study of combustion of natural gas 
combustion in cylindrical shaped turbulent diffusion flame 
furnace is presented in this paper. The objective of this study is 
to predict the NOx content in natural gas combustor based on 
species transport model. The mathematical models used for the 
combustion of natural gas and air includes standard k-t model 
for turbulent flow and species transport model. The effect of 
mass flow rate of air and methane, stoichiometric ratio and mass 
fraction of air on the flame shape is study in this paper. The 
result is observed in this study of species transport model on the 
basis of furnace performance with different equivalence ratio. 
For this study of combustion process, mesh of working 
combustion model developed on commercial software program 
FLUENT was imported to simulate the data. The cylindrical 
shaped combustor is used to burn the methane -air using eddy 
dissipation/finite rate eddy dissipation model. The model used 
for one step combustion reaction of methane-air assuming 
complete combustion of the fuel to carbon dioxide and water. 

Index Terms —thermal, prompt, natural gas combustor, 

I. Introduction 

Combustion is one of the most important processes widely 
used in many industrial applications. The study of combustion 
has reached their advanced technology in the past few 
decades. For the prediction of NOx formation three different 
mechanism of NO formation is used, that is thermal, prompt 
and fuel nitrogen conversion. NO is generally formed by three 
chemical routes. Zeldovich or thermal mechanism that is 
depend on the temperature where rate formation of NO is 
exponentially dependent on flame temperature below 1850 K. 
prompt or Fenimore mechanism generally having complex 
chemistry and it increases with fuel rich combustion as 
compared to lean combustion. Our aim of this study is to 
analyze the mixing of chemical species in a cylindrical shaped 
turbulent diffusion flame furnace. Many researchers worked 
in this field via experimentally as well as computational way 
to analyze the different combustion process. One of them is 
cited here: K.M.Pandey, D.H.Das and B.Acharya worked on 
2D model CFD analysis of combustion process. In their 
research, they concluded the effect of variation of specific 
heat on temperature in gaseous combustion with Fluent 
Software. 

Mean flow Equation: 
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Where t^is the mean velocity in the direction X ■ ,p the fluid 
density, pi the viscosity,the prandtl number of the chemical 
species,i?jthe rate of generation of the chemical species per 
unit volume ,K the thermal conductivity, Q-react the 
volumetric rate of heat generation , tj r(3d the volumetric rate 
of heat transfer with radiation and u t u } and U t h are the 
unknown Reynolds stresses and heat fluxes, respectively. 

II. Numerical Model 

A cylindrical shaped combustor mesh is imported for the 
simulation of methane-air. The fluent 15.0 solves the 
governing equation of continuity, momentum, species and 
energy in the gas phase. The one step global mechanism used 
for methane air reaction is given as below: 

CO 2 ~h 2.H n 0 T- / .52 

The mechanism shown above for methane -air reaction 
contains 5 species and 1 reaction is employed in the present 
study. 

Based on the following governing equations, 2D simulations 
are performed. The finite volume method are used for the 
discretization of governing equations and SIMPLE 
algorithms are used is used for pressure and velocity coupling. 

The equations are solved implicitly with a 2D pressure based 
solver using an under relaxation method. 


Fuel Inlet Velocity 

0.4 m/s 

Air Inlet Velocity 

lOOm/s 

Equivalence Ratio 

0.80 

Excess air percent 

25% 

Inlet temperature 

300 K 
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III. RESULT AND CONCLUSION: 


Table. 1: Values observed after simulation of the 

DATA 


Exit temperature 

602.75 K 

Exit velocity 

0.814m/s 

Mass fraction of pollutant 
NO assuming only 

thermal NOx 

0.002 

Mass fraction of 

pollutant NO assuming 
only prompt NOx 

0.00016 

Mass fraction of pollutant 
NOx assuming thermal 
and prompt NOx 

0.0023 


As per literature the NOx concentration is increases with 
increasing equivalence ratio at first and then decreases. In the 
fig . 1 the contour of mass fraction of NO ppm is shown, which 
is calculated by assuming thermal and prompt NOx together. 
Similarly, the contour of the mass fraction of NO ppm 
assuming thermal and prompt NOx separately are shown in 
fig.2 and fig.3. NO concentration is very high where 
temperature is very high and having nitrogen and oxygen 
concentration are available there. Mass weighted average exit 
NO mass fraction is lower without prompt NOx as compared 
to thermal and prompt NOx together. So we have also 
examined that in this case flame is lean and prompt NO 
production is low. 
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Contours of Mass fraction of Pollutant no Dec 19, 2014 

ANSYS Fluent 15.0 (axl. pbns, spe, ske) 

Fig.l Mass fraction of NO ppm assuming thermal and prompt 
NOx together 
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Dec 19.2014 

ANSYS Fluent 15.0 (axl. pbns, spe. ske) 


Fig.3 Mass fraction of NO ppm assuming prompt NOx 
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Fig.2 Mass fraction of NO ppm assuming thermal NOx 
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